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The magnetic response of strongly localized electrons to a time-dependent vector potential is 
considered. The orbital magnetic moment of the system, away from steady-state conditions, is 
obtained. The expression involves the tunneling and phonon-assisted hopping currents between 
localized states. The frequency and temperature dependence of the orbital magnetization is analyzed 
as function of the admittances connecting localized levels. It is shown that quantum interference 
of the localized wave functions contributes to the moment a term which follows adiabatically the 
time-dependent perturbation. 
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I. INTRODUCTION 

The explanation of the orbital magnetic moment at 
thermal equilibrium in terms of "ring^purrents" has been 
known for more than half a centuryEI. The most stud- 
ied configuration is that of a small mesoscopic structure 
penetrated by a constant magnetic flux, in which the 
electrons move ballistically or diffusivelyo. Recently, the 
problem of the ring currents and the associated mag- 
netic moment have been addressed in the case of strong 
disordeiu. In this regime the electronic states are local- 
ized and their overlaps are small. It has been shown that 
since the coupling of the electrons to the thermal bath 
enables resonant processegfl, there appears a current flow- 
ing in opposite sense to the one due to tunneling, which 
exists in the absence of the electron-phonon interaction. 
This counter|-current has been related to the dc Hall effect 
in insulatorsO. 

The study offthe orbital magnetic response to a time- 
dependent fluxEHZl is of interest by itself and also because it 
is relevant to experimentaa deigned to detect persistent 
currents in mesoscopic ringsS. Nonequilibrium magne- 
tization may he used to measure specific electronic re- 
laxation ratesEEI. The magnetic moment in response to 
a time-dependent flux has been discussed so far in the 
context of weak disorder. 

In this paper we investigate the response of strongly 
localized electrons to a time-dependent flux, which sets 
up a magnetic moment in the system. We obtain the 
nonequilibrium orbital magnetization and study it when 
the electrons are coupled to the thermal bath. 

The calculation is based on the Kubo approach that 
yields, for strongly localized electrons, the current be- 
tween two localized states in response to the time- 
dependent flux. The proper definition of the orbital mag- 
netization away from steady-state conditions, in terms of 
these currents, is a somewhat delicate point. We devote 
the next section to a discussion of this issue. In Sec. Ill 
we obtain the current that flows between two localized 



states. Section IV contains the analysis of the character- 
istic frequencies that dominate the nonequilibrium mag- 
netization of an isolated group of three localized states. 
(This is the smallest cluster which can have an orbital 
moment.) Section V contains concluding comments. 

The formalism we develop allows for the calculation of 
the ac conductance of strongly localized electrons coupled 
to a thermal bath. We derive in the Appendix quantum 
interference corrections to the ac hopping conductance 
and show that they yield nonlocal contributions to the 
'bond' admittances. 



II. NONEQUILIBRIUM MAGNETIZATION OF 
LOCALIZED ELECTRONS 

Our aim is to obtain an expression for the orbital 
magnetic moment of localized electrons, at nonequilib- 
rium conditions. We will show that the moment can 
be given in terms of the currents between the localized 
states. When these currents flow in response to a con- 
stant magnetic fleld, i.e. at equilibrium, this is a straight- 
forward task. In that case the equilibrium value of the 
orbital magnetization, Meq, is just a manifestation of 
the Aharonov-Bohm effect and can be obtained from the 
thermodynamic potentiaB. The result reads 



M,. 



(1) 



in which lij is the current between the state localized at 
Ri and the state localized at Rj, and {£j) indicates a sum 
over pairs. At equilibrium, current conservation implies 
that the 'bond' currents lij satisfy Iij = 0, which in 
turn yields that M^q is independent of the choice of the 
coordinate origin. At nonequilibrium one has 



Pj 



(2) 
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where pj is the time derivative of the electron occupa- 
tion on the jth locahzed state. The contribution of the 
dielectric currents resulting from the nonequilibrium oc- 
cupations has to be included in the orbital magnetic mo- 
ment. Eoi' an homogeneous continuous system, this is 
achievedliil by introducing a single vector function that 
encompasses both the (spatial-dependent) charge and the 
current densities. In the present case, the analogous 
quantity is 



where TV is the number of localized states, 
because the generalized currents Jij satisfy 



(3) 



This is so 



(4) 



which follows from the requirement of charge conserva- 
tion Pj = 0). Note that the requirement (^) corre- 
sponds to a current density with zero divergence in the 
case of a homogeneous continuous system. As a result, 
the nonequilibrium orbital magnetic moment, M, is given 
by 



M 



(5) 



<ej> 



and is also independent of the choice of the coordinate 
origin. It takes a particularly convenient form when the 
relative coordinates ff 



re 



Re — i?, R 



(6) 



are introduced {R is the coordinate of the center of grav- 
ity). Then, using Eqs. (0) and (0) we obtain 



M 



(7) 



In this expression the radius-vectors Re appearing in (|i|) 
are replaced by the relative coordinates re- It is seen 
that for the choice (||) the displacement currents disap- 
pear from the expression for M. 



III. THE LINEAR RESPONSE EXPRESSION FOR 
THE CURRENT 

The properties of strongly localized electrons interact- 
ing with a phonon bath are most conveniently-discussed 
using the Hamiltonian introduced by Holsteinllj. In that 
description the electron site density is coupled to the 
ion displacements. (The term 'site' denotes a localized 
state.) The Holstein Hamiltonian, in the presence of a 
vector potential, can be written in the form 



n = Y^ ^tc\ce + '^g^g-*? + H JejQeje'^''c\cj, (8) 



in which Qej results from the electron-phonon interaction 



hj = exp 



(9) 



Here ce and are the annihilation and creation oper- 
ators of an electron in the state localized on site I and 
te are the on-site energies, assumed to be randomly dis- 
tributed (energies are measured from the Fermi level). 
The phonon operators are denoted and 6t, and uj^ is 
the phonon energy. In Eq. (^), Jij is the overlap of two 
electronic wave functions localized on sites j and and 
is the matrix clement for the electron-phonon interaction 
on site The strong localization regime is characterized 
by \Jij\ ^ ~ The transfer of an electron between 
two sites, as described by the last term of 7i, can be ac- 
complished by tunneling, or can be accompanied by the 
exchange of one or more phonons with the thermal bath 
(as is evident by expanding Qej in terms of v^J ). 

In this description the vector potential applied to the 
system appears as a phase factor, expi(/)£j, in which (pej is 
the time-dependent flux (in units of the flux quantum 0o) 
contained in the area [Re x Rj\/2, where Ri is the radius 
vector to the site I. The effect of this flux is treated as a 
perturbation. Accordingly, we re-write the Hamiltonian 
(P) in the form 



Ti = Tic + Ti 



ph 



Til + Hp 



(10) 



in which 



Tic = ^ ^ec\ce, Hph = ^ ^9*^91 
e q 

Hi — ^ JejQijclcj, 

refer to the unperturbed system, and 

Tipor = ^ JejQe]i(l)e]{t)c\cj. 
to 



(11) 



(12) 



The current flowing in response to the timfedependent 
flux can be obtained by using the expressiora 



ije = 2eQJejQeje"'"^clcj 



(13) 



for the current operator. Then linear response theory 
yields that the net current between sites j and £ is 



I At) = 



dti 



npcr{h),2'iJejQej{t)4it)cj{t) 



-cl>ej{t){2mejQej{t)4{t)cj{t)). 



(14) 
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In this expression the time dependences of the opera- 
tors, as well as the thermal averages (indicated by the 
angular brackets) are determined by the Hamiltonian 
Ho + Hph + Hi, Eqs. dnj). 

We now derive an explicit expression for the current, 
valid in the strong localization regime. As this regime 
is characterized by the small parameter \Jje\/\f-j — q|, 
we expand Ij£ in powers of the wave- function overlaps. 
Clearly, [cf. Eqs. ( p^ ) and (p^)] the expansion starts at 
order J^. Quantum interference effects appear firstu in 
order J^, which allows for the interference of the direct 
amplitude j —>■ ( with the indirect one, via a third site 



k. Such processes lead to the appearance of "ring cur- 
rents" for a group of at least three, sites in the case of 
a constant (time-independent) fluxH. In this paper we 
concentrate on effects resulting from the second order in 
the overlaps. However, to demonstrate that quantum in- 
terference effects of strongly localized electrons can be 
studied using the Kubo approach, we present in the Ap- 
pendix the derivation of the bond current to third order 
in the overlaps. 

To second order in the wave- function overlaps the bond 
current is given by 



e\Jje\^(bej{t) / dtJn.il - ni)2^e''''^'-''\Qjeiti)Qij{t)) - n,(l - n,)23e'^«(*-*i) (g,,(i)Q,,(ii))) , (15) 



where we have used (pjg — —(pij. (The superscript (2) 
indicates a result to order J^.) Here nj denotes the equi- 
librium electron occupation of site j, Uj = {e^'^j + 1)^^ 
and egj = ei — ej. 

The time-dependent correlators of the phonon opera- 
tors Qj£ are calculated with respect to the free phonon 
Hamiltonian, Tiph. One finds 



{Qij{t)Qji{ti)) 



9(t) = E - 



[1 + 7Vg-)e-'''^f* + iV,-e'"«-* 



(16) 



where iV^- = (e^'^« — 1) ^ is the Bose function. Here 
we have taken vtv^_^ ^ 5ji\v^\'^, which is basedEl on 
?;J. = v^expiiq - Re). We now insert Eq. ([l^ ) into Eq. 
( |l5|) and put 



where (jjjg is the amplitude of the ac flux 
follows that 



~(2) 

with the amphtude I^^ given by 

X Lj-(1 - n^)25e"«*+^»(-*)-29(o) 

~ne{l - nj)23e"«*+29(*)~2s(o)l ^^g^ 

~(2) 

One notes that I^^ {uj) vanishes as the frequency tends 
to zero. In other words, to order the current is driven 



(17) 
It then 

(18) 



I);>{uj)^e(j3ji\Jji\^ j dt(l~e- 



by the time dependence of the flux. This is not the case 
in higher orders in the overlaps. We show in the Ap- 
pendix that due to quantum interference of the electron 
wave functions (which appear first to order J^), the bond 
currents and consequently the orbital magnetic moment 
have a part that follows adiabatically the external flux. 

(2) 

One can now define an effective admittance zj^ for the 
j£ bond 



(2), 



(20) 



For weak electron-phonon coupling, the real part of the 
admittance becomes 



(2) 



(21) 



with 



1 

2^ 



(1 



+ {e^'^ - 1 



dte'^'^'+'^'>\l + g{t)) 
dte''-'^'-'^^\l + git))], 



lijii^)- 



(22) 



(It is assumed that the localized levels are not acciden- 
tally degenerate.) It is seen that two types of transitions 
between the levels £ and j contribute to , resonant 
transitions and non-resonant, phonon-assisted ones that 
exist due to the coupling with the thermal bath. At low 
frequencies such that (3uj 1 we have 

SRzj^-' = e^l3\Jj£fnj{l - nii)iT S{eji + uj) + 5{e.ji - uj) 

+ Gj,, (23) 
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with Gji being the dc bond conductance- At low tern- iujCjS^j — zgj (ioj—H 
peratures such that /3|e| ^ 1, it becomesEJ £ V 2c 



[re X rj\ 



ejl + kfl + lej-d 



(24) 



Here v is the interaction matrix element for the electron- 
phonon coupling, and P is the phonon density of states. 
The resonant transitions, described by the first term in 
( p3| ) require that the probability to find a resonating pair 
of sites will not vanish. In a bulk system this mecha- 
nism may therefore be expected to contribute. The full 
quantum-mechanical analysis of an isolated pair of sites, 
valid for arbitrary frequencies, is given e.g. in Ref. n4. 

(2) ... . — 

The imaginary part of z]^^ is associated with the dielec- 
tric response of the j£ bond, and is out-of-phase with the 
real term. In the absence of the coupling to the thermal 
bath it becomes 



2 rij — n£ uj 



(25) 



The coupling with the phonons introduces additional 
contributions to the dielectric response which are ignored 
here for simplicity. 



IV. ANALYSIS OF THE MAGNETIC MOMENT 

In exploiting the Kubo formula, one should interpret 
the driving force as an effective one, which arises from 
the externally applied field and the fields induced in 
the system. This means that the time derivative of the 
flux shouldJncorporate the changes in the local chemical 
potentialsEj. Using the relative coordinates introduced 
in ^ one thus has 

— c -> 

iujcfije — iuj—H ■ [fj x fi] 

-iLu — [H X R] ■ {fj - fi) + Snj - S^ii, (26) 

in which Sfij is the change of the chemical potential at 
site j. The first term here results from the time depen- 
dence of the magnetic field H. In fact, {e/2c)H ■ [rj x ri] 
is just the magnetic flux (divided by the flux quantum 
^o) acquired from the magnetic field between sites j and 
£. The second term can be interpreted as arising from an 
electric field, E = iuj[H x R]/2c. Accordingly, the last 
three terms in ( |2^ ) represent the local electro-chemical 
potential difference. 

To determine the local chemical potentials, and con- 
sequently the bond currents, we need to solve the set of 
equations ((2|), where 



Pj — iLue(3nj(l — nj)5p,j 



(27) 



and lij is given by Eq. (^0|). (We omit the superscript 
(2) for brevity.) Thus we have 



iloyW X R] ■ in - Tj) + 5pL( - 5^jij ) , (28) 



where 



Co =e2/37ij(l-nj) 



(29) 



is the capacitance connected to the jih site. As the cal- 
culation is restricted to terms linear in the magnetic field 
H the second term on the right hand side of (^ can be 
omitted since it vanishes upon averaging over the site lo- 
cations. It follows that the local chemical potentials and 
consequently the currents are determined by the mag- 
netic fluxes Lpej 



(30) 



The magnetic moment, Eq. (|^, is given in terms of 
the bond currents. To find those, we use Eqs. (0) to 
write [c/. Eq. 



ituCj 



(31) 



We then obtain 



lt-2 ( ■ 

e V 



iujCe^ 



(32) 



This can be rearranged to yield 




V/, 

2^1 



Ztlf.2 T 

iLoCi^ e 



(33) 



Equation (^3|) represents an array oi C x C equations for 
the bond currents, where C is the number of bonds in the 
system. We shall not attempt to present the formal gen- 
eral solution for this array. Instead, we will obtain the 
limiting behaviors for the frequency dependence of Af , 
and exemplify them by considering the smallest possible 
cluster of sites, a triangle, which possesses a magnetic 
moment. (For a two-site cluster the orbital magnetiza- 
tion (0) vanishes because [ri x ^2] = by construction.) 

Inspection of Eqs. (^3|) reveals the appearance of the 
characteristic frequency cJq , given by 



0^ 



C. 



(34) 
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When the external frequency lo is much higher than all 
luq 's, the local chemical potentials are ineffective in de- 
termining the bond currents. In that case 



hj ^ — fijZij. 



(35) 



On the other hand, when w is less than all w^^'s, the lo- 
cal chemical potentials dominate the bond currents. In 
that situation [cf. Eq. the current in the j£ bond 

"feels" the fluxes on all other bonds, and the solution of 
Eqs. (|3^) is much more complicate. For example, for a 
three-site cluster we find 



_ _ tOJif f 1 1 ^ \ 

Jl2 — ^23 — ^31 — I 1 1 I 

e Vzi2 ^23 Z3i/ 

where ip is the total flux enclosed by the triangle 



^ = ¥^12 + </'23 + 



(36) 



(37) 



That is, ia;(y9/c is the e.m.f. applied to the cluster. Recall- 
ing that for a triangle r i x r 2 = ^2 x ra = f 3 x ri — 2S/3, 
where S is the (vectorial) area of the three-site group, 
we find that the magnetic moment can be written in the 
form 



-^triangle — ~ ^ ' 



(38) 



where Zj^ stands for the largest admittance of the three in 
the high frequency case, and for the smallest one at very 
low frequencies. Note, however, that for a triangle, the 
numerical coefficient in the high frequency limit includes 
a geometrical factor ((^12 /St/?) which diminishes the ratio 
M/w. 

Let us now focus on the frequency dependence of the 
admittance itself. We consider the cases where oj <C \ejt\ 
and Lo 3> One can then ignore the resonant contri- 

bution to 3?Zj£. For uj <^ \eji\ we have 



OJ 



Zje = Gjtil + i—) 



with bji^ given by 



e 2 



(39) 



(40) 



In the other extreme limit, where u ^ Icj^l, the imag- 
inary part of the admittance becomes inversely propor- 
tional to UJ. When the electron-phonon coupling is weak, 
Zji in this regime is dominated by the imaginary part, 
leading to 



low frequencies, such that ui < \€ji\ and oj < ujI^, the ad- 
mittance is mainly given by the bond dc conductance, so 
that 



-^triangle ~ ^ ^ 3 



(42) 



The magnetic moment decays exponentially as T 0. 
As the frequency increases such that lo > w^^, then 
iGjgUj/ujl^ and 



triangle ' o 



^HS''G,i/uji'. 



(43) 



The leading frequency dependence is modified as com- 
pared to the very low lu regime, but the moment is still 
exponentially decaying with the decrease of the temper- 
ature [see Eqs. (|4|) and (^o|)]. Finally, in the high fre- 
quency limit such that lu ^ \£je\, the moment becomes 
almost independent of the temperature and the frequency 



triangle 



(44) 



For frequencies in the intermediate regime the behavior 
is rather complicate, and depends upon the detailed rela- 
tionships between the frequency, temperature, \eji\, and 
the characteristic frequencies uj^^ and uj{^ . 



V. CONCLUDING REMARKS 

The orbital magnetization of localized electrons in re- 
sponse to the temporal variation of a vector potential has 
been studied. We have obtained an expression for the 
magnetic moment in terms of the tunneling and phonon- 
assisted hopping currents between localized states. These 
currents have been derived from the Kubo formula. 

The mere existence of an orbital magnetic moment in 
the strongly localized system requires a group of at least 
three states. We have found that the magnetization of 
this small cluster is proportional to the relevant domi- 
nant admittance between sites (the smallest one for small 
frequencies, uj ^ uj^^, and the largest one in the oppo- 
site case). Consequently, Mtriangio is proportional to | Jp, 
where J is the overlap integral of two wave functions. The 
equilibrium magnetic moment of a triangleB, on the other 
hand, is of order J^. This is because Meq results from 
quantum interference processes among the three sites. 

The magnetization of a bulk system involves summa- 
tion over all bonds [see Eq. (^]. Its full calculation re- 
quires the knowledge of the distribution function of the 
site energies and the localized site locations. This may 
lead to a temperature dependence different from the ex- 
ponential one characterizing the small cluster. Consider 



Zji 



iZji/uj, Zji ^ e'^\Jji\'^{ne - nj)/eej- (41) for example the case of nearest-neighbor hopping 



Collecting the above results, we obtain the following 
pattern for the behavior of Mtriangio as function of the 
temperature and the external frequency. Starting at very 



at 

low enough frequencies. Then the typical energies in- 
volved are within a region of width ~ /?~^ around the 
Fermi level. Consequently, for a smooth distribution of 
the site energies, within an energy band of width W, the 
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number of triangles contributing to the magnetic moment 
is oc {PW)~^. The typical spatial separation between 
sites is ^ ^Lp, with ^ being the localization length and 
Lp (X In 1/3 Jo I . (Jo is the bare overlap integral between 
two localized levels.) It follows that the area of a typi- 
cal triangle is of order (^L^)^. Collecting these estimates 
and using Eq. (^2|) we obtain that the bulk magnetic 
moment per unit volume is 

Mb^ik ^ tu;{-rHv'V{P~^)(3{^)^^L0 fe- (45) 
c pW 

Here rii is the concentration of localized levels. The bulk 
moment is proportional to the frequency and depends 
on the temperature T as T'^V{T)\n^ {Jq/T). Modifica- 
tions, due e.g. to the temperature dependence of the 
relevant area in the case of variable-range hopping, are 
possible. Finally we rempjk that in this work electronic 
correlations were ignoredllj. These have to be taken into 
account in the calculation of the bond currents. We hope 
to pursue this issue in a future article. 
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APPENDIX A: QUANTUM INTERFERENCE 
EFFECTS IN THE CURRENT 



As is mentioned above, quantum interference between 
the direct amplitude j £ and the indirect one, via a 
third site k, appears first in order J^. It is therefore nec- 
essary to derive the corrections to the current to third 

(3) 

order in the overlaps. We denote this quantity by I]i ■ 



Returning to Eq. (14) we find 



dti / dt2 



t .t 
dti / dt2 



-(l)ij{t){ I dti I dt2 (^i(i2)^i(ii)23fiJ£jg£j(0c](i)cj(i) + 2»J<;jQ^j(i)4(t)cj(i)^i(ii)^i(i2 

— OO 'J —OQ 

dti f dt2(niiti)2nji,Qi,{t)4it)cjit)ni{t2))\. (Ai) 



ac amplitude of I'^f^ , I^f^ , the expression 



The calculation of the terms appearing in Eq. (Al) is straightforward, though rather cumbersome. It yields for the 



OO /'OC 



-^j?('^)=2e^ / dTi / dT2<FM,kji-Ti,T2) 



+ <Pkee 



+ (Pike~ 



nkUjUl - Fjk,tk{-Tl,T2) 



x\cbe,+^,ke-'^'' +cl>kee-''^^- 



rijUkUi 



+Fjl^ki{-Tl,-T2) 



+4>kie 



-iuj{Ti+T2) 



niTljUk + Fij^kjiTl,T2) 



rijUkUi - Fej^ikiTi,T2) 



+ c^.ke-''^^' 



+ (l>kie-"^'' +<pjke 



niUkUj 



(A2) 



where rij = 1 — nj, and 



^>j,fc,(Tl,T2) = K{jjfeJfc,J,,e"«^i+^^'=^^^+s(^i+^^)+s(^i)+3(^^)-33(0)y 



(A3) 



Contrary to the behavior of lj'f\uj), — 0) is finite. Setting w = in Eq. (A2) one finds that each fc— term in 

the sum is proportional to 



r(3)r,., - 



(A4) 
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that is, to the total flux enclosed in the triangle £jk. The zero frequency part of the current is just the persistent 
current flowing in response to a|-jConstant magnetic flux. Its properties in the strong localization regime have been 
studied in great detail elsewhereEl. Denoting 



(A5) 



we derive from (A2) I-^ (lj) — Ipc for weak electron-phonon interaction, and in the small lj limit. The result is 



lf{u^)-I,c 



JjkJkeJe^ 



1 1 i« 

! jfJuinj^gj 



^ e^i ^kj 
2iLU 



Pjk 



1 fUk -Ui Hi - rij 



kl ~ + 

^fej ^ ^Ik 



Pjk 



ki £j£ 

1 fT-k — rij rij — Hi 



f-lk 



njUkJjk 



(A6) 



A remarkable observation is that l\f^{u)), which is part 
of the current in the bond, is driven by the magnetic 
phases on the bonds and jfc, but not by the phase on 
the bond j£ itself. This is a manifestation of interference 
effects. The interference affects other electronic proper- 
ties of the system. One example is the conductivity. The 
formalism presented here can be used to find that quan- 
tity. In that case the electric field E is related to the 
fluxes by 



lujq 



eE ■ R 



■jg- 



(A7) 



Then the ring currents Ipc vanish by virtue of (A41. Us- 
ing Eq. (A7) one finds that Eqs. (pO[), and (A6) give 
the conductivity of the j^-bond. In particular, Eqs. ( pOj ) 
and ( p^ ) yield the well-known result for the bond con- 
ductivity to second order in the overlaps, while Eq. ( [A^ ) 
gives the rwulocal corrections to it, which result from 
interferenceH'tl. Namely, the current in the j-^-bond is af- 
fected by the voltage drops on the k£ and the jk bonds. 
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